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1 Introduction to BFKL 



Hadronic processes at high-energy colliders often involve more than one energy 
scale. As a consequence, calculations in perturbative Quantum Chromodynamics 
(QCD) can involve large logarithms of the ratio of these scales, which must be re- 
summed to obtain a reliable prediction. Two processes where this might be necessary 
are Deep Inelastic Scattering (DIS) at small x and hadronic dijet production at large 
rapidity intervals Ay. In DIS the logarithm that appears is ln(l/x), with x ~ Q 2 /s 
the squared ratio of the momentum transfer to the photon-hadron center-of-mass en- 
ergy. In large-rapidity dijet production the large logarithm is Ay ~ ln(s/|t|), with 
s the squared parton center-of-mass energy and |t| of the order of the squared jet 
transverse energy. In both of these cases the large logarithms, which arise at each 
order in the coupling constant a s , can be resummed by means of the Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) equation [[[]]. 

The most familiar prediction of the BFKL resummation in the leading logarithmic 
(LL) approximation is the power-law rise in the partonic cross section as a function 
of the energy: 

a ss e AAy S3 s A . (1) 

The quantity (1 + A) is often referred to as the BFKL Pomeron intercept, where at 
LL 

A = a s 41n2 (2) 

with a s = N c a s /ir and N c = 3 the number of colors. Similarly, in DIS the structure 
functions are predicted to rise as x~ A at small x. One of the goals in BFKL physics 
has been to observe this power-law rise as a direct indication of the importance of the 
resummation. In this talk I will review the status of BFKL physics, both in its theo- 
retical development at next-to-leading logarithm (NLL) and in its phenomenological 
application to experiment. 



2 BFKL at LL 

I begin by giving a simple physical picture of the BFKL resummation at lead- 
ing logarithm (LL), where the large logarithm is taken to be the rapidity interval 
between two widely-separated partonic jets. Although this in no way can be consid- 
ered a derivation of the BFKL equation, it is useful to show what assumptions go 
into the resummation and to show how the factors of a s Ay arise at each order and 
exponentiate. 

The starting point is the factorization of the partonic cross section at large rapidity 
separation: 

d^ a %p b± = vM ^ K A v>P^>r^ v b ( P l L ) . (3) 
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Figure 1: A schematic picture of the cross section for producing particles at large rapidity 
separation. 

A physical interpretation of this factorization is represented in Fig. |l|. The process 
consists of two distinct scatterings, which occur at widely-separated rapidities, y a and 
y b , and small transverse momenta \p a ±\ ~ \Pb±\- The impact factors V(p 2 ± ) depend 
only on the transverse momentum exchanged and the specific partons involved in 
each scattering, but not on anything else that happens in the event. The function 
f(Ay,p a ±,p bl ) connects the two scatterings by accounting for the emission of gluons 
in the rapidity interval. This function is universal and provides the exponentiation of 
the logarithms. 

We can see this factorization directly in the Born cross section for gluon-gluon 
scattering at y a y b , shown diagrammatically in Fig. ||(a): 
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where the central factor is the O(a s ) contribution to the function /, and we see that 
the leading-order gluon impact factor is 
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The real 0(a).) correction to / can be obtained by considering the emission of three 
gluons, strongly ordered in rapidity y a ^> y\ y b , shown diagrammatically in 
Fig. 0(b): 
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From this formula, we easily see where the large logarithm comes from. It arises from 
the integral over the rapidity y\ of the intermediate gluon, resulting in a factor of 
a s Ay. 

We can now generalize the form of these real corrections to an arbitrary number 
of emitted gluons: 

• At each order, another gluon is inserted in the ladder with a weight given by 

a s d 2 k i± 

TT kf ± 

• The emitted gluons conserve transverse momentum, as enforced by a delta- 
function factor 

• The intermediate gluons are integrated over the ordered rapidity-intervals, 

Vb < Vl < Vl < ■ ■ ■ < Vn < Va , 

producing an overall factor of 

n\ 

The real n-gluon contribution can be obtained directly from the tree-level (n + 2)- 
gluon cross-section by assuming that all transverse momenta are comparable in size, 
while terms suppressed by 0(e~^ Vi ~ y ^) are neglected. 

Of course, the real contributions by themselves are not infrared-safe, because they 
diverge when the k± of any of the intermediate gluons vanishes. The cure for this is the 
inclusion of virtual corrections in the large Ay limit. The 0{a 1 s ) virtual correction 
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in gg scattering comes from the diagram in Fig. |2|(c). Similar contributions must 
come in at each order to regularize the l/k\ infrared singularities. Including these 
corrections, we can understand the following form of the BFKL equation 



f(Ay,p aX ,Pbx) = t:<5 (2) \p a ± + Pb±) 
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where the kernel K is an integral operator acting on a function <p(p a ±) by 

d 2 k 
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At LL, each operation of the kernel inserts one more real or virtual gluon in the 
ladder. The first term in (^) corresponds to a real gluon, while the second term is the 
virtual correction needed to regularize the soft singularity. The equation (0) is just 
the solution of the BFKL equation obtained by iteration. 

A more compact form of the BFKL solution at LL is obtained by finding the 
eigenf unctions and eigenvalues of the kernel (M). This yields 



f(Ay,p a± ,p bl _) 
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a s x(n, 7 ) = a s 2^(1) - V(f + 7) " V»(f + 1 " 7) 



(9) 



(10) 



gives the eigenvalues of the LL BFKL kernel (H), and ip is the logarithmic derivative of 
the gamma function. For very large Ay, the integral over 7 in (0) can be performed 
in the saddle-point approximation. The n = term dominates, and one obtains 
the exponential rise in the cross section (|I|), displayed in the introduction with A = 
a a x{0A) = 4a s ln2. 



3 BFKL at NLL 

At LL each application of the kernel gives a contribution of 0(a s Ay). At NLL 
one also includes terms of 0(a 2 s Ay). That is, we can reinterpret the kernel in (0) as a 
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power series K = a s K^ + a 2 s K^ + . . ., where a s K^ l > is the LL kernel, given in (j8|), 
and a 2 s K^ includes the NLL corrections. There are three types of contributions at 
NLL, which are shown schematically, in the context of gg scattering, in Fig. fj. They 
consist of: (a) the emission of two gluons nearby in rapidity, (b) the virtual correction 
to the emission of one gluon, widely separated in rapidity, (c) the subleading purely- 
virtual corrections. These three types of contributions took many years and many 
papers to sort out the technical details^, with the final NLL kernel obtained in 1998 by 
Fadin and Lipatov |J. Although the full kernel has not been checked in a completely 
independent manner, many of the pieces of the calculation have received independent 
confirmation. Two particularly significant checks are the calculation of the virtual 
correction to the gluon emission at large rapidity separation || , and the compilation 
of the three NLL terms into a single kernel with the cancellation of all collinear and 
soft singularities |>j. 

The final result of this calculation is usually presented by applying the NLL kernel 
to the LL eigenfunctions, with azimuthal averaging, yielding 



K 



NLL 



(PL) 7 " 1 {a^MT)}^) 7 ' 1 

a s (X>X (1) (7) [l - a a (fj,)b ln(plj > 2 ) 

+ *.(aOV%)) (p^) 7 " 1 (11) 



where a 8 X (7) is the LL eigenvalue for n = given in (|10l), &o — 11/12 — n//(6iV c ), 
and we have explicitly included the dependence on the MS renormalization scale \i. 
The NLL correction has been separated into two terms. The first term depends on 
the scale p a ± and is associated with the running of the coupling in the LL kernel: 
— > a s (p a ±). The second term, a^x^Kl), is independent of scale and contains 
the remainder of the NLL corrections [[J. 

After completion of the NLL corrections to the BFKL kernel, several issues quickly 
became apparent. Roughly speaking, they can be separated into issues associated 

1 A list of references can be found in ref. H, but with no guarantee of completeness. 
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Figure 4: Leading BFKL eigenvalue A = a s x(k) at LL and NLL. 



with the running coupling term and issues associated with the scale-invariant term. 
In this talk I will concentrate on the scale-invariant term. For analyses of the running 
coupling issues, see Refs. ||[7],[8||| . 

The first indication of problems with BFKL at NLL was seen immediately by 
Fadin and Lipatov. The corrections to the leading eigenvalue are large and negative! 
If we ignore the effects of running coupling, we obtain 



«*x(i) 



2.77a, - 18.34a" 



(12) 



for three active flavors. This function is plotted in Fig. f|. At the not-unreasonable 
value of a s = 0.16 the NLL corrections exactly cancel the LL term, while for larger 
values of a s the eigenvalue becomes negative. Naively, this would indicate that the 
BFKL Pomeron intercept also becomes negative, leading to a cross section that de- 
creases, rather than increases, as a power of the energy. 

Unfortunately, things get even worse. The standard BFKL power-law scaling 
of the partonic cross section ([1]) relies on the saddle-point evaluation of the NLL 
generalization of the BFKL solution @. Upon closer analysis Ross |JT0] showed that 



the NLL eigenvalue function x{l) no longer has a maximum at 7 = ~, but has 
a minimum with two maxima occuring symmetrically on either side of this point^. 
Performing a higher-order expansion of x(l)i Ross found a smaller correction to the 
BFKL Pomeron intercept. However, the cross section he obtained was not positive 
definite. It contained oscillations as one varied p a ± and pbj_- This led Levin [§] to 
declare that NLL BFKL has a serious pathology. 

2 The standard procedure in these analyses is to modify the LL eigenfunctions used in eq. ([ll]) in 
order to make the eigenvalues manifestly symmetric under 7 — > 1 — 7, following ref. S. 
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One might wonder whether the approximate evaluation of the integral performed 
by Ross is adequate at this stage. Perhaps an exact evaluation is necessary. However, 
negative cross sections have also arisen when the resummed small- a; anomalous di- 
mensions, obtained from the NLL BFKL solution, were used to study DIS scattering 



at small-a; |[ll| , p!2| . In any event the NLL corrections to the BFKL solution are large, 



leading one to question the stability and applicability of the BFKL resummation 
procedure in general. 



4 Understanding the large NLL Corrections. 

When any perturbation expansion has large corrections at higher orders, the nat- 
ural thing to do is to try to reorganize the series so that it converges more rapidly. 
In this talk I will briefly discuss three different approaches to this reorganization. 
The first proposal by Brodsky et al. |13|] uses the freedom to choose the renor- 



malization scheme. The NLL eigenvalue equation (|i~2[) is written in the MS scheme. 
Brodsky et al. argued that a non-abelian physical scheme should be more natural 



for the BFKL resummation. Then they used the BLM procedure |14| to find the 
optimal scale for the QCD coupling. In this case the BLM procedure dictates a large 
scale, thereby reducing the effective a s (and the LL prediction) and also reducing the 
coefficient of the NLL a 2 s term. This approach predicts a BFKL Pomeron intercept 
of A ~ 0.17 for a s = 0.2. In addition it yields a very weak dependence on the gluon 
virtuality p^ ± and leads to an approximate conformal invariance. 



The motivation for the second proposal (first suggested in |15j and |16 1 ) can be 
seen from the discussion of the physics of BFKL at LL in section The approximation 
used in deriving the LL contribution at each order in the BFKL ladder was to neglect 
terms of 0(e _ ' yi_yi+1 ') in the QCD matrix elements, which is valid when the emitted 
gluons are all widely separated in rapidity. However, the gluon rapidity yi is then 
integrated all the way up to Thus, the errors in the matrix elements are largest 
when yi ~ This suggests that one enforce a condition y i+1 — > A, so that 

the gluons are required to be widely separated, and the kinematic approximations 
are good. The arbitrary parameter A is assumed to be much smaller than the total 
rapidity interval. The excluded region is re-introduced at NLL, such that the change 
in the cross section due to shifting A is always next-to-next-to-leading logarithm 
(NNLL). In this way, the dependence on A can be regarded as an estimate of the 
uncertainty due to NNLL corrections (similar to the role of the renormalization scale 
/z in the MS scheme). 

Fig. [| shows the dependence on A of the leading eigenvalue and its second deriva- 
tive at LL and NLL for a s = 0.15 in this modified BFKL theory. Note that the 
corrections to a s x{\) are not large for A > 2 and have weak dependence on A for 
large A. Also, the point 7 = ~ is a maximum for this coupling as long as A > 2.2. 
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Figure 5: Dependence of A = ctsxi^) an d B = ~2 as ^" (2) 011 ^ ^ or as = 0-15, from 
Ref. §. 



Thus, the BFKL resummation is stable for large enough A. This modification of the 
BFKL resummation predicts a somewhat larger value of the BFKL pomeron intercept 
than the previous proposal. However, the implications of a large value of A for the 
phenomenological use of BFKL is open to interpretation. 



The third and perhaps most ambitious proposal |17|,|9| to control the large NLL 
corrections is to systematically include the largest collinearly-enhanced, but still sub- 
leading corrections into LL BFKL. The most important of these are energy-scale 
corrections 0. To understand the origin of these corrections, note that in our dis- 
cussion of BFKL at LL in section ^| we chose to work with the symmetric rapidity 
Ay = In s/(p a ±pb±) as the large logarithm to resum. However, we could equally well 
have chosen y + = \nx^/x^ = Ins jp\ or y~ = lnx^/x^ = lns/p^, where xf is the 
momentum fraction along the positive or negative light-cone for the emitted gluon i. 
These choices are all equivalent at LL because the transverse momenta are treated as 
comparable in size; however, at NLL they are inequivalent. A change in the logarithm 
produces a change in the NLL kernel and can introduce double transverse logarithms 
of the form a s In (Pa±/pl±) the resummation. 

Motivated by DGLAP-type resummation one finds that the appropriate choice 
is to resum y + when pf ± 3> p 2 aL and y~ when p 2 aL 3> pl ± . The effect of these 
changes of the BFKL resummation variable was studied in refs. [||] and ||, and 
the corresponding terms in the NLL eigenvalue x il) were isolated and resummed 
in ref. ||17|| . Additional collinearly-enhanced terms due to the effects of the running 
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Figure 6: The BFKL Pomeron intercept (here labeled u) at LL, NLL and with the 
collinearly-enhanced resummation included at NLL, from Ref. ||. The curve oj s corresponds 
to the coefficient for the exponential rise in Ay, which is what we have focused on in this 
talk. The curve cj c corresponds to the power growth of the small- x splitting functions. 
Although identical in LL BFKL, these two coefficients differ at NLL. 



coupling and the non-singular part of the splitting functions have also been considered 



and resummed in Ref. ||. A nice discussion of these ideas can be found in ref. [19| . 

A similar approach, advocated in Refs. |[20|| , is to use the "duality" relations 
between the BFKL ln(l/a;) and the DGLAP ln(Q 2 ) resummations to incorporate 
the dominant collinear effects into BFKL at small x. Adding contributions obtained 
from the known LO and NLO DGLAP anomalous dimensions, one gets a "double 
leading" expansion for the BFKL function xil)-, which is better behaved and more 
stable in the collinear region near 7 = 0. 

Results of the collinearly-enhanced resummation from Ref. || are shown in Fig. || 
where the leading eigenvalue u> s = a s x(|) is plotted as a function of a s . As in the 
other cases the eigenvalue is found to be positive after resummation, yielding a value 
of about u s = 0.27 for a s = 0.2. In addition the characteristic eigenvalue function of 
7 is stable for values of a s of interest. 

The physical implications of the energy-scale dependence can be seen by further 
investigating the relation between resummation in y^ and Ay. When p% ± S> p^±, 
the resummation in y + requires the ordering > x£ . Translating back into the 
symmetric variable, this implies Ay > ln(p 6J _/p a _|_). Similarly, when p 2 a± 3> plj_, the 
resummation in y~ requires the ordering x^ > x~, implying Ay > ha.(p a ±/Pb±)- These 
constraints hold for any two successively emitted gluons. Therefore, the incorporation 
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of these collinear effects corresponds to imposing a p±- dependent cut, yi + i — yi > 
| ln(pij_/pj + ij_)|, on the separation in rapidity between the neighboring gluons. Since 
this cut is very similar to the rapidity veto of proposal two, it is understandable that 
when both the collinear resummation and the rapidity veto are included, as studied 
in Ref. the dependence on the parameter A was significantly reduced, even for 
small A. 



5 Phenomenology of BFKL 

Although the theoretical studies of BFKL have been focused on its behavior at 
NLL, the level of most phenomenological studies is still at LL. To treat a process 
consistently at NLL, one also must incorporate the 0(a s ) corrections to the impact 
factors. Although these corrections are known for some processes, at least at the 
amplitude level, they have yet to be incorporated into a consistent calculation suitable 
for phenomenological studies. 

In this section I will discuss several probes of BFKL physics in hadron-hadron, 
lepton-hadron, and 7*7* collisions. I will only consider processes for which the rele- 
vant transverse scales p a _L,Pb±. on both sides of the BFKL ladder can be considered 
perturbative. In particular, I will not consider the inclusive -^(x, Q 2 ) in DIS, because 
it is necessarily dependent on non-perturbative inputs. 

5.1 BFKL probes at the Tevatron 

One of the most thoroughly studied searches for BFKL physics has been by the D0 
collaboration in pp collisions at the Fermilab Tevatron. As in all BFKL experimental 
studies, the basic idea is to analyze the events in a configuration which most closely 
approximates that used in the BFKL resummation. Jets, with transverse momentum 
above some -Ej_ m i n , are tagged and ordered in rapidity. Then one defines observables 
as a function of Ay = y a — yb, where y a and y^ are the rapidities of the most forward 
and backward jets, respectively. 

The most natural BFKL signal would be the power-law growth in the partonic 
cross section with Ay, as in Eq. (H). However, at fixed center-of-mass energy this 
growth is swamped by the effects of steeply falling parton distribution functions 
(PDFs) which are relevant when far forward or backward jets are produced. Thus, 
the first observable to be considered was the decorrelation in azimuthal angle between 
the two tagged jets as a function of Ay. Physically, this effect is easy to understand. 
In the Born approximation, only two jets are produced, and by momentum conser- 
vation they must be back-to-back. However, as the rapidity interval increases, there 
is more room for additional jets, and the tagged jets become decorrelated. This can 
be seen directly in the BFKL solution (|9]). At small Ay, all terms in the Fourier 
series in (f> = (p a — fa — ir are approximately equal, producing a delta-function in (ft 
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which forces the two jets to be back to back. As Ay increases, the higher order terms 
become smaller and smaller compared to the leading n = term, so the jets become 
completely decorrelated p2| . The simplest observable to display this effect is the 



moment (cos0), which goes to 1 if the jet azimuthal angles are completely correlated 
and goes to if they are completely decorrelated |23|] . 



Even before the D0 analysis was completed, however, it was realized that there 
was a serious problem in using LL BFKL for phenomenological analyses at hadron 
colliders: the BFKL resummation includes the contribution of energetically disfavored 
or disallowed configurations in its predictions. In principle these configurations are 
subleading, but in practice they are very important fl24| . (In fact these effects could 



be considered to be a foreshadowing of the large corrections to BFKL at NLL.) We 
can understand this effect by considering the Feynman x-values used in the PDFs. 
The exact values are given by conservation of light-cone momentum along the beam 
axis and can be written 

x a = -j= (p a ±e Va + p b± e Vb + J2 h±e y ^j 

x b = -j=( Pa ±e- ya +p b ±e~ yb + l>^ r% ) > ( 13 ) 



where the sum is over all partons produced in the event. In the "naive" LL BFKL 
one only keeps the leading contributions, 



Pa±e Va 



H = — r~ ■ ( 14 ) 
V s 

That is, one convolutes the analytic LL BFKL solution (|9|) with the impact factors, 
using the PDFs evaluated at x Q ab . However, the true x a ^ and are always larger than 
x®b, an d the energy-momentum constraints x a> t> < 1 are not enforced in the naive 
BFKL calculation. If the PDFs vary strongly with x, this can greatly overestimate 
the contributions from multi-jet events. 

With the analytic LL BFKL solution (|9|) the phase space of the intermediate gluons 
has already been integrated over, so there is no choice but to use the leading x°'s (0) 



in the PDFs. However, in a BFKL Monte Carlo solution [ p5|j26[| one generates the 
gluon ladder directly as in Eq. ([7|). Thus, one has information on all the produced 
partons, and one can enforce energy conservation on the solution by using the exact 
x's ( |l"3"D in the PDFs. This greatly improves the reliability of the BFKL prediction. 

In Fig. we show the D0 azimuthal decorrelation data from Ref. [B7J compared 
with the naive LL BFKL and a Monte Carlo BFKL calculation with energy con- 
servation included. The data shows (cos</>) as a function of Ay at a/s = 1.8 TeV, 
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Figure 7: Azimuthal angle decorrelation as a function of rapidity interval. The D0 data is 
from the lowest curve is "naive" BFKL, the two upper curves are BFKL Monte Carlo 
predictions. The correlated systematic error is due to jet energy scale uncertainty. 



where the tagging jets are required to have E± > 20 GeV. The naive BFKL severely 
overestimates the rate of decorrelation with Ay, due to the overweighting of energeti- 
cally disfavored or disallowed configurations. The two BFKL Monte Carlo curves are 
produced by the Monte Carlo in Ref. [25). Both use the exact x's in the PDFs, but 



use different approximations to relate the partonic to the hadronic cross sections in 
the high energy limit. The difference between the two approximations is subleading 
in the high energy limit. The BFKL Monte Carlo with energy conservation certainly 
works much better than the naive BFKL, but the subleading uncertainties are still 
sizeable, as displayed by the difference between the two Monte Carlo predictions. 

An argument against the decorrelation measurement as a signal of BFKL is that 
it probes the region around <fi ~ 0, which is also most sensitive to Sudakov logarithms. 
It would be nicer to probe directly the rise in partonic cross section with the partonic 
energy as in Eq. (|l|). This became possible when the Tevatron collider was run at the 
lower energy of 630 GeV, allowing a comparison of dijet production at two different 
center-of-mass energies. By binning the events in the partonic x a ,x b values rather 
than in Ay, the dependence on the PDFs should cancel in the ratio 

= 1800GeV) 



cr 



= 630GeV) 



R 



This is the original proposal of Mueller and Navelet 
point approximation, one obtains a prediction of 

e A(A yi -Ay 2 ) 

%FKL = 



Ayi/Ay 2 



(15) 

Using the asymptotic saddle 



(16) 



where Ay { « ln(si/p a ±p b ±) « ln(si/E] 



This ratio was measured by the D0 



collaboration, in several x bins, with the following cuts: E 



±mm 



> 20 GeV, \y\ < 3, 
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Ay > 2, and 400 < Q 2 = E a± E b± < 1000 GeV 2 . Using the formula (pD, a BFKL 
Pomeron intercept of A = 0.65 ± 0.07 was extracted . 

This measurement is noteworthy in that it is probably the only current mea- 
surement which shows a rise in the cross section that is larger than the LL BFKL 
prediction. Using a s (20 GeV) = 0.17 in Eq. one obtains a LL prediction of 
A = 0.45, which is almost 3 standard deviations below the extracted value. However, 



as discussed in Ref. one must be careful in interpreting this measurement. First, 
the extraction of the BFKL Pomeron intercept from (pl|) assumes that the asymp- 
totic BFKL expression is valid, and the experimental cuts and precise definition of 
the x's do not significantly affect the asymptotics. In particular the cut on Q 2 was 
seen to slow the approach to asymptotics, resulting in a smaller predicted value for 
the ratio R. The inclusion of energy conservation via a BFKL Monte Carlo further 
reduced the predicted ratio. Finally, it was shown that the use of equal E± cuts on 
both of the tagging jets introduces the same large Sudakov logarithms that plague 
the decorrelation measurement. Thus, it seems unlikely that the large ratio found by 
the D0 collaboration can be attributed to perturbative BFKL. 

5.2 BFKL probes at HERA 

It is also possible to look for the BFKL rise in the cross section in deep inelastic 



scattering (DIS) by tagging a forward jet [ 3~I] , j32]| . Referring to the high energy factor- 



ization picture of Fig. [I], the DIS setup consists of the scattering of an off-shell photon 
with virtuality Q 2 and Bjorken x b j = Q 2 /s ep on the left, the scattering of a forward 
jet of momentum fraction Xj e t = E- ]ct / E p and transverse momentum p± on the right, 
connected by the BFKL ladder of gluon emissions in the middle. If Q 2 ~ p 2 ± ^> A.q CD , 
then the BFKL evolution is perturbative. The large logarithm that is resummed is 

ln(XjetMj)- 

This DIS setup has several advantages over the pp setup, due to its asymmetric 
nature. Note that in the high energy limit, the PDF of the proton is evaluated 
at Xjet- Thus, with Xj e t fixed, one can vary the rapidity interval at a single collider 
energy by varying x b j, without any change in the PDF. In addition, this suggests that 
the energy conservation effects mentioned above may be less important, or at least 
not so strongly dependent on x b j. Finally, the resummation in ln(xj et / 'scy), rather 
than in the jet rapidity intervals, is more natural here, and perhaps more stable 
theoretically, since it is corresponds to the standard DGLAP evolution variable when 
Q 2 > p\. The only major disadvantage is that the virtual photon impact factor is 
more complicated theoretically than the gluon or quark impact factors. Indeed, it 
has not been calculated completely at NLO. 

Both the HI |33| and ZEUS [[34]] collaborations have measured this forward jet 



cross section. The main experimental cuts on the forward jet itself are x- jet > 0.035, 
£jL je t > 3.5 and 5 GeV for HI, x ]et > 0.036, E ±iet > 5 GeV for ZEUS, and 0.5 < 
E±jet/Q 2 < 2 for both. The HI data is displayed in Fig. |[ In Figs. |](a) and (c) the 
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Figure 8: Forward jet cross-sections as a function of Bjorken-x from HI [^] for two p± 
cuts of 3.5 and 5 GeV. The data in (a) and (c) are compared with Monte-Carlo model 
predictions, including hadronization. The data in (b) and (d) are compared with partonic 
NLO 0(a 2 ') and LL BFKL calculations. 



data is compared against several hadron-level Monte Carlo shower models attached 
to lowest order QCD matrix elements. The two that agree the best with the data 
are ARIADNE |5[ and RAPGAP f3§. ARIADNE is based on the colour dipole 
model for gluon radiation which, like BFKL, lacks ordering in k±. Gluon radiation 
in RAPGAP is based on DGLAP evolution, but this model also includes a resolved 
photon contribution to the basic QCD production mechanism. In Figs. |8](b) and 
(d) the HI data is compared against a LL BFKL prediction |57] and a NLO QCD 
0(a 2 s ) prediction [pq| , both at the parton levels. The NLO calculation significantly 
underestimates the data at small Xbj, whereas the BFKL calculation overestimates it. 
This is not unreasonable, given that NLL corrections to BFKL are expected to reduce 
the rise at small x^, and that the kinematic cuts could not be included exactly in 
the calculation. The ZEUS data [34], shown in Fig. P, similarly is far above a NLO 
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Figure 9: Forward jet cross-sections as a function of Bjorken-j; from ZEUS (34|] for p± > 5 
GeV. The data are compared with several Monte-Carlo model predictions at the hadron 
level. 



QCD calculation [pS|| , but below the LL BFKL expectations. A fit [4(| to the data 



from both experiments using the LL BFKL cross section yielded an effective BFKL 
pomeron intercept corresponding to A = 0.43 ± 0.025(stat) ± 0.025(sys), compared 
to the LL prediction from Eq. (@) of A = 0.75 for a s = 0.28 at Q 2 = 10 GeV 2 . 

Recently, several forward jet calculations with different approaches have shown 
good agreement with the data. One calculation [^T| is based on LL BFKL, but 
modified by a consistency condition containing effects similar to the dominant NLL 
energy-scale effects discussed in section |[ A second calculation |42| is with a hadron- 
level Monte Carlo generator, based on the CCFM evolution E3|, which is designed to 



agree with both DGLAP and BFKL in their respective regimes of reliability. Since 
both of these calculations can be considered LL BFKL, with some dominant subiead- 
ing corrections included, this looks promising for BFKL. A third approach which 
also fits the data is a NLO calculation that includes a resolved photon contribution, 
similar in spirit to the RAPGAP Monte Carlo. Interestingly, it appears that the 
success of this approach relies not on the evolution in Q 2 allowed by the inclusion of 
the photon PDF, but in the fact that it effectively approximates one term higher in 
a s , via the NLO resolved piece. This is not incompatible with BFKL, since the new 
0(a^) contribution also includes the first gluon emission in the BFKL ladder. It is 
an interesting question to ask how the approximations in this picture of forward jet 
production mesh with those in the BFKL picture. 
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Figure 10: Two-photon cross-sections, ct^*^*, after subtraction of the LO 'y*^/* — * QQ 
contribution at y/a ~ 189-202 GeV and 4 < Q{ 2 < 40 GeV 2 , from Ref. The top figure 
compares the data with LO BFKL in the saddle point approximation with A = ap — 1 = 0.53 
and with the one-gluon exchange diagram. The lower figure shows fits of the saddle-point 
BFKL expression to the data, allowing either ap or the normalization K to vary. 



5.3 BFKL probes in 7*7* 



Another standard BFKL measurement |45[ is to observe the total 7*7* hadronic 
cross section as a function of ^/s 77 . This can be extracted from the process e + e~ — > 
e + e~+ hadrons by tagging on the forward and backward electrons. In this case the 
two independent scatterings of Fig. [I] are the off-shell photons of momenta Q\ and Q\ 
that break up into color dipoles (qq pairs at leading order), which are then connected 

by the BFKL ladder. The large logarithm here is Y ~ ln(s 77 / ^ Q1Q2) , i n direct 
analogy to the hadron-hadron case. 

This experiment has the advantage that, for large enough Q\ 2 , there are no non- 
perturbative PDF inputs, so that in principle it is the cleanest probe, theoretically. 
However, from a purely calculational point of view, this process may be more com- 
plicated since it involves two off-shell photon impact factors. In particular, the cross- 
channel gluon on which the BFKL ladder is built does not even appear until NNLO 
in a standard perturbative QCD calculation.. 



The 7*7* hadronic cross section has been measured by the L3 p6|j47|] and the 
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OPAL |48] collaborations at LEP. In Fig. [H| we display the most recent preliminary 



L3 data [ i7j . The data clearly shows a rise with Y as expected by BFKL, but much 
less steep than LO BFKL. However, the data is above both LO and NLO QCD 
predictions In addition the data is also above the prediction from the single- 

gluon exchange contribution (evaluated using the average value of (Qf 2 ) = 15 GeV 2 ). 
In the lower figure the data is fitted to the asymptotic BFKL prediction 
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(17) 



with either the overall normalization or the BFKL intercept A = ap — 1 left as a free 
parameter. The preliminary result in the latter case gives A = 0.36 ± 0.02, which is 



more in line with NLO expectations. The OPAL preliminary measurements j|8[ give 
qualtitatively similar results, but with less statistical significance. 



6 Summary and Conclusions. 

In this talk I have presented some recent results in BFKL physics, both in the 
theory and in its phenomenological applications. On the theoretical side, the focus 
has been on the large NLL corrections. At this time, it seems fairly safe to say that 
the original catastrophe of falling, or even negative, cross sections has been averted. 
By understanding the origin of these effects in the collinear behavior of the gluons, 
one can reorganize the resummation in order to move the dominant corrections back 
into the LL theory. Then the NLL prediction for the BFKL intercept is stable and 
slightly smaller than the standard LL prediction. 

On the phenomenological side, the results look suggestive, especially in the DIS 
and 7*7* data. The cross sections are significantly above the state-of-the-art NLO 
QCD calculations, as they should be for a BFKL-enhanced observable. They are also 
in the range one might expect from a NLL BFKL calculation. However, the lesson 
learned from the analysis of the hadron-hadron experiments is that one must be 
very careful to consider how experimental cuts and kinematic effects, such as energy 
conservation, will affect these predictions. Although the asymmetric configuration of 
the DIS forward jet experiments and the lack of nonperturbative PDF-dependence in 
the 7*7* (assuming Q 2 is large enough) may make these observables less susceptible to 
large subleading effects, a thorough phenomenological analysis is certainly warranted. 

With the completion and understanding of the NLL corrections to the BFKL 
kernel in hand, the next phase is to bring the phenomenological analyses up to the 
same NLL level. So far the major emphasis has been on the NLL BFKL intercept, 
but to make a full NLL prediction, with a reliable normalization, one also needs to 
combine this with NLO impact factors. In particular the NLO impact factor for 
the off-shell photon is crucial for the both the DIS and 7*7* analyses. In addition 
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the analyses must be performed in such a way as to treat the kinematics and cuts 
as accurately as possible. Promising steps in this direction are Refs. fn)JI7|, which 



incorporate the largest NLL corrections into BFKL via a consistency condition, or via 
a CCFM Monte Carlo, respectively Another useful exercise would be to incorporate 



the largest NLL corrections to the BFKL ladder into the BFKL Monte Carlos [25,26 



and to modify them for use in the other experimental environments. This would 
be helpful for gauging the sensitivity to subleading kinematic effects (which arise at 
least as much from the impact factors as from the actual BFKL ladder). However, 
the greatest progress would come with a full NLL calculation with full NLO impact 
factors included. There is much work left to be done. 
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